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Abstract 

Heavy hadrons are analyzed in a random and dilute gas of instantons. We 
derive the instanton-induced interactions between heavy and light quarks at 
next to leading order in the heavy quark mass and in the planar approxima- 
tion, and discuss their effects on the hadronic spectrum. The role of these 
interactions in the formation of exotic hadrons is also discussed. 
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I. INTRODUCTION 



Hadrons with one or many heavy quarks exhibit a new type of symmetry: invariance 
under spin-flip of the heavy quark Q. This invariance can be used to organize the hadronic 
structure and properties of heavy-light systems. A number of relations follow both in the 
spectrum and among form factors of heavy hadrons when the mass of the heavy quark is 
taken to infinity M. 

The interplay between light and heavy degrees of freedom in heavy-light hadrons can be 
clarified in the heavy-quark limit by combining chiral symmetry with heavy quark symmetry 
||. The basic observation is to note that the hard part in the heavy quark field is kinematical 
and factorizable. The remaining part is soft and constrained by chiral dynamics. Hence, 
the soft physics in heavy-light systems can be analyzed in a way similar to the light-light 
systems. A qualitative understanding of this part can be achieved by using QCD inspired 
models. 

In this paper we discuss the effects of a random gas of instantons and antiinstantons 
on mesons and baryons containing one or several heavy quarks. We analyze the correlation 
functions of various hadrons with one or many heavy quarks in inverse powers of the heavy 
quark mass mq using a succession of Foldy-Wouthuysen transformations prior to radiative 
corrections. In section 2, we give a brief summary of the salient properties of some typical 
heavy-light systems. In section 3, we show how the heavy meson correlator may be system- 
atically analyzed in inverse powers of the heavy quark mass, in the planar approximation. 
Recoil and magnetic corrections to both the heavy quark propagator and the heavy-meson 
correlator are evaluated at next to leading order in the heavy quark mass. These results are 
quantified in the form of effective interactions between heavy and light constituent quarks. 
In section 4, we estimate the effects of the induced effective interactions and heavy meson 
and baryon spectra. The results are in overall agreement with the constituent quark model 
and heavy solitons. We briefly discuss similar effects in exotic hadronic configurations. Our 
conclusions are summarized in section 5. 
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II. GENERALITIES 



Throughout, heavy hadrons will be understood as mesons or baryons with at least one 
heavy quark. Although we will be interested in taking the heavy quark mass to infinity 
and then relaxing it, we will in fact specifically have in mind for a heavy quark, the bottom 
quark b with a mass of 4.7 — 5.3 GeV, the charmed quark c with a mass of 1.3 — 1.7 GeV, 
and to some extent the strange quark s with a mass of 100 — 300 MeV. 

Heavy mesons Qq may be organized in their ground state into multiplets with I(J P ) = 
|(0~, In the heavy quark limit, the multiplets are invariant representations of the heavy 
quark symmetry group (essentially left spin rotation). Empirically (K,K*) = (493,892) 
MeV, (D,D*) = (1869,2010) MeV and (B,B*) = (5278,5324) MeV. D* decay is mainly 
through Dtx. 

The heavy baryons with one heavy quark will be of the type Qqq. The conventions being 
Aq = 0| + , Eg = l| + and Eg = l| + . The ones with two heavy quarks will be of the 
type QQq. Few heavy baryons have already been observed. Organizing them in multiplets 
invariant under heavy quark symmetry, we have A s , (E s , E*) = 1116,(1195,1385) MeV, 
A c , (E C ,E*) = 2284, (2455, 2530 ±5 ±5 @Q) MeV, and A 6 , (E 6 , EJ) = 5641, (..., ...) MeV. The 
dots refer to yet to be measured masses. Other measured heavy (charmed) baryons include 
(S c , H*) = (2468, 2642.8 ± 2.2 §) MeV and fi c = 2704 MeV. Heavy baryons with more than 
one heavy quark have not been found yet. 

The basic principles at work in a heavy light system are best illustrated using a simple 
bag model description. If we were to insert a heavy source in a spherical cavity of radius 
R, then the total energy can be organized using the bare heavy quark mass mq following 
E = rriQ + Eqitlq + EiTTIq 1 + .... The contribution E refers to the energy of the light quarks 
present in the cavity and is standard |J. The contribution E\ corresponds to 



To be confirmed by other experiments. 
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m Q 2m Q 2m Q 

where the first term is the recoil of the heavy quark, and the second term is the magnetic 
interaction between the average magnetic field induced by the light quark at the center of the 
bag, and the magnetic moment of the heavy quark. While schematic, ([!]) captures the essence 
of the 1/m.Q corrections in heavy quark physics. Using standard bag model parameters ||, 
we have for charmed mesons (recoil, spin) ~ (400, 20) MeV, while for bottom mesons (recoil, 
spin) ~ (100, 5) MeV 0. 

III. HEAVY HADRONS IN AN INSTANTON GAS 

In what follows, we will try to understand the origins of Eq and E\ from a microscopic 
description of the QCD vacuum using a random gas of instantons and antiinstantons. 

A. Heavy Quark Expansion 

Consider the correlation function of a heavy-light meson in the QCD vacuum. In 
Minkowski space, it reads 

C±(x,x') = (T (qT ± ij(x)^T ± q(x')) > (2) 

with r± = (1,7) x (1 ± 7°)/2 x (1, T) a non-relativistic source with arbitrary flavor. Here 
q(x) refers to the light quark, and i/j(x) to the heavy quark with bare mass rriQ. For mq 
much larger than the typical scale of the problem Aqcd, one may use Aqcd/^q expansion 
to analyze (|2]). We perform this expansion using a Foldy-Wouthuysen transformation of the 
heavy quark field ||. 

where V = d — iA and with g (the gauge coupling) set to one. The first transformation 
rescales the momenta, the second eliminates the odd parts, and the third removes the mass 
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term. The successive transformations in (2) are vector-like, unitary and gauge-covariant. In 
terms of (0) the QCD part of the action for the heavy field ip becomes 

~ Qii°V°Q + QOiQ + 0{m Q 2 ) (4) 

where operator 0\ is defined as 

( 5) 

2m Q 2m Q 

with B l = — ie y "*[V*, V J "]. Equation (Q) has the expected FW form to order nig 1 . In the 
first term in (H) we recognize the recoil, in the second we recognize spin effect on the heavy 
quark in external field . Under (^|) the heavy meson source shifts 

As a result, the correlator (0) takes the generic form 

iy- V 



QQ, 



r?|o> (7) 



2m Q 

Mixing between the particle and the antiparticle content of the correlator ([?]) drops out in 
the nonrelativistic limit 0. Thus, in Euclidean space 

C±(x, x') ~ - e ^(^') (0|Tr (r ± Sq(x, x') Y ± S(x', x)) |0) (8) 

where S is the propagator of the light quark, and Sq is the heavy quark propagator, 

S Q ~ + SooOiSoo + 0(m Q 2 ) (9) 

with 5*00 = 74/^4 being the free part. This construction can be carried out to arbitrary 
orders in 1/toq HQ, given that the heavy quark expansion is not upset by renormalization 



B. Heavy Quark Propagator 

The heavy quark propagator ([5|) may be analyzed in a random instanton gas. In the 
planar approximation, the infinitely heavy quark propagator satisfies the integral equation 
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s^ 1 = s; 1 -Y,((4-,}-s oo y 1 ) (io) 
1,1 

where ^4 = 7 4 A 4 and = ry 4 <9 4 . The sum is over all instantons and anti-instantons and 
the averaging is over the position zj and the SU(N C ) color orientation Uj, with 

A I {x)=U I A{x-z I)P )U\ (11) 

Generically 



2 \VJ 'J J V / 2V 4 iV c _ 

where Tr^ stands for a trace in color space. For a low instanton density n = N/V4 ~ 1 
fm~ 4 , we can iterate (|10"D in powers of n. The result is 

s^ 1 = s- 1 + m ( e + — e! + 0(^-2-)) + ( n2 ) ( 13 ) 

where the diluteness factor is given by the dimensionless combination np 4 ~ 1CT 3 . Substi- 
tuting (|HJ) into (|TDD , we obtain 



6 = I d 4 zj Tr c [ S; l (—^— - S*)S; L + I I ) (14) 



■1.1 



1= I^TiclS: 1 !^ - Si )o l (^—-SAS; 1 + I^I) cir,) 



and 



where Ci = — V 2 /2 — cr • 5/2. Both O and 0i are r— dependent. To proceed further, we 
note that in coordinate space, the heavy quark propagator in the one instanton background 
reads 



< x 



-J_|0 >= 5{x) 6{t) l±J- p e */ T *M*.-*i) (i 6 ) 
% V4 / 2 



with x s = (s,x). Inserting (|TB|) into (]TJ],]n^) and using the one-instanton configuration 

a;(x) = +% vXv (1 _ _L_) (17) 



yield for large times 

< x-oc|(e , Qi)\x +00 >~ 87rp 3 (-4I , +Ii) (18) 

with 



rio,ii, 



f 00 x 2 dx {cos\f x /2) , ({d^osf x ] 2 -[d^mf x ] 2 + { P A^) 2 C os2f x )} (19) 

J 



where f x = ir\x\/y/l + x 2 . There is no contribution to 0i from the spin part —a ■ B. From 
([L8D it follows that the instanton induced shift in the heavy quark mass is AqMq ~ 70 MeV 
pT0|| from Go and AxMq ~ 16 MeV from Qi for a c— quark with m c = 1350 MeV. The recoil 
effect AiMq is an order of magnitude down compared to the naive bag estimate ([I]). This 
can be understood by noting that in the presence of instantons, the energy of a heavy quark 
can be rewritten schematically as 

E = 32tt x no 4 x (l/p+ + ..\ (20) 
V m Q J 

where the factors follow from (|13|) and (fl8|) . (p0[) is the analog of (fj). For the instanton 
parameters used, and a charmed quark, (|20|) yields 

E ~ 327r x 10~ 3 x (^600 + - 600 + ~ (^60 + 30 + MeV (21) 

which shows that the zeroth order shift in the mass is about 60 MeV, while the recoil effect 
is about 30 MeV, as expected. 

Similar arguments can be used for the light quark propagator. The result is S^ 1 ~ 
Sq 1 + i^fnYjiyx) with an average light quark mass shift AM, ~ 420MeV P,|TT . 



C. Heavy Quark Correlator 

The correlator (fj), written generically as C ~ (S 1 ®^), receives contributions from both 
planar and non-planar graphs and is usually hard to analyze in the random gas approxima- 
tion exactly. In the planar approximation, things simplify. After resummation, the inverse 
correlator @j reads 
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c- 1 ~ s- 1 ®sl- 1 -J2((s-4 J 1 )- 1 ® OS 



oo -7 1 4 / 



(22) 



The upper script T is short for transpose. Here the spin-flavor-color indices are left uncon- 
tracted and the space-time indices are omitted. Standard perturbation techniques in the 
massless sector gives 

1 



4-i ~ s 











s~ l -4i 


\®n) 



(23) 



where |$ n ) is the normalized eigenstate of the Dirac operator in a one-instanton background, 



(i0 - 4 /) $„(x - z{) = \ n $ n (x - Zl ) 
It follows that the 't Hooft zero mode $>o(x — %i) |0 

1 p a; • 7 



(24) 



$o(ac) 



7T (l 2 + p 2 ) 3 / 2 v^ 2 

is dominant in (|23|). Thus, to leading order we obtain 



(25) 



V" 1 / 



c 1 ~ S" 1 <g> 5, 



1 o cT,-l 



n / dSj Tr c ([ L ]/ <g> [ H ]j + J -»• / 



(26) 



with 



S"o ) Sn 1 



|$o)($o| 
i-v/n£ / 

t-jL — - s.) s; 1 + — s; 1 f-^— - 5,) 0! 

K «7 V 4 ,/ / rriQ \J7 4 V 4 ,j / 



(27) 



z 7 4 V 



where S = ($ |S|$ ) ~ (240 MeV)" 1 . 



D. Effective Interactions 



The inverse correlator (|26| ) allows for an immediate translation to effective interactions. 
In the long wavelength limit, the instanton size is small, and a local interaction between the 
effective fields Q and q can be derived much like the 't Hooft interaction between the light 
effective fields q JTSyrjJ. From the Bethe-Salpeter equation associated to (fffiD, we read the 
vertex 



S 



r%{x,y, x',y') = -inN c j d 4 Zj J dU 1 (t/f (x\[Lj })\x')U^ b ® U c k (y\[ Hi }t\y')U^ + I - /) 

(28) 

where the color matrices have been explicitly displayed. This vertex function gives rise to 
an effective action Si 

T a b c d (x, y, x', y') = 1 — 1 (29) 

8q»(x) 5qt(x>) 8Q*(y) 8Q\{y>) 

which is essentially non-local. In the long-wavelength approximation the kernel (^8]) factor- 

izes into two independent kernels as x — > x' and y — > y'. The corresponding effective action 

reads 



Si = —inN c J d A zi J dUi J d A x c^(x) Ui (x\ Li \x) llj q(x) 



d 4 y Q^y)Ui(y\Hi\y)UjQ(y) 



+ 1^1 (30) 



£! Q = n(- 1 -^)(^-) fiQtil^Qiqtq+^iQti^A'QiqtA^ (31) 



and yields the effective interaction in Euclidean space (leading order in 1/N C 

? 3 I Q \ / 4ttV 
N c /VVnSoA^ 2 ~* "* M ' 4"^ 2 

For the detailed construction of the above lagrangian we refer to the Appendix and our 
previous paper ( ||). The first bracket in ( pl|) arises from the heavy quark part and the 
second bracket from the light quark part. Wick-rotating to Minkowski space gives 

which is to be compared with the 't Hooft vertex for two light flavors q = (u, d) 

/ AM 2 \ 

£ qq = ( ) (det q^qi + det q L q fi ) (33) 



Interaction ( p2|) is dominated by the Coulomb-like second term and has a proper heavy 
quark spin symmetry. The recoil effect renormalizes the strength of the interaction through 
AMq = AqMq + AiMq ~ 86MeV. The spin part gives rise to a chromomagnetic interaction 

AM a AM" pin 1 _1 + 7° 
™ = 9 Q_ i Qi+J_ X a huq qAV^q (34) 

qQ 2nN r 4^2 V ' 



with 

AMn in = — n 167rp 3 / x 2 dx ^f* ~ 3 MeV (35) 
w prriq J (1 + x 2 ) 2 

for a c— quark. As a Coulomb-like term in fl32|), it has a smooth Af? limit for the large N c 
and is attractive in the spin zero, color-singlet channel. 

Similar arguments may be applied to the heavy mesons QQ as well. To order 1 / rriQ and 
in the planar approximation, the effective interaction among the heavy quarks is given by 

/ AMn AMn \ / 1+7° —1 + 7° 1 — 1+7° —1+7° \ 

The recoil effects appear to first order in I/tuq and renormalize AMq. The spin effects are 
of second order in I/tuq, and result in 

A£3? = ( AKryT ) iQi±2! A vrq Q^Avrq (37) 

For heavy baryons of the type qqQ we have 

/ AMqAM 2 \ /_1 + 7 ° 
£<?<?<? = 2n2Ar2 J 2 Q ( det qLq/J + det q/? q L ) + 

\ Q^^Vq (detq L A a q fi + detq^A^ ) ) (38) 

and to second order in l/rnq 

/AMT n AM 2 \ 1 _1 + 7 ° / \ 

4*3 = " ( n 2 N 2 ) I Q — g-^-AV^Q detq L A < Vq R + detq R A V^q L J (39) 

The phenomenological implications of these interactions on heavy-light spectra will be dis- 
cussed next. 



IV. HEAVY HADRON SPECTRA 

The contributions of the various instanton interactions derived above to the heavy hadron 
spectra, will be discussed using a variational approach. For Mesons, the generic Hamiltonian 
is 
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where M is the reduced mass of the heavy-light system, with m q = AM q ~ 420 MeV 
and rriQ = m c + AMq ~ (1350 + 86) MeV. The harmonic potential provides for a simple 
mechanism of confinement. The instanton-induced interaction derived from Eqs. ( $2\ 
38|) will be treated as a perturbation. The trial wavefunction is 

1>(X) = (^) e"^ 2 (41) 

where \ = (f q — tq). Minimizing the expectation value of (f5]) in (f41~D with respect to a 
yields a = \Muj with the confining energy E a = |u; as expected. Since the size r = 
of the ground state is a function of the reduced mass M, we fix our parameters by the size 
of the heavy-light system r q q = 0.6 fm. Then the size of the the heavy-heavy system is 
tqq ~ 0.4 fm, and the confining energy is about £ a ~ 250MeV for both of them. 

A. Mesons 

For heavy mesons the relevant two-body instanton-induced interactions are for Qq 

r(2) _ f AMq AM q " 



For QQ 



tt(2) 
n QQ 



V 2nN r 



(AMqAMc 



V 2nN c 

The induced spin-interaction in the Qq configuration is given by 

' AMn in AM a \ 1 



") (i + ^^-rg) (42) 



^ + \^Q 2 )^(r Ql -r Q2 ) (43) 



n 2 ' s 



^-a q -a Q \ a q .\ Q 5\r q -r Q ) (44) 



qQ " V 2nN c 

We recall that AMq 111 is down by one power of l/irtq. Using fl42|) and the trial wavefunction 



(flT|) we have 

, f2 v f AMnAM„\ . ,^.„ N c (AM AM g \ ( 1 \ 3 

n 



and similarly for QEj). Thus ( H$ ) 183 MeV and ( H { q ] q ) ~ -103 MeV. These 

numbers should be compared respectively to —140 MeV and —70 MeV, as quoted in || 
using qualitative arguments. The spin corrections are of order ( H^q ) ~ —42 MeV, a result 



that is consistent with the constituent Quark Model estimate of ~ —27 MeV [15]. The spin 



induced interaction in heavy systems such as charmonium is tiny, (Hqq) ~ —0.7 MeV. 

To evaluate the spectrum in heavy-light and heavy-heavy systems, we use the mass 
formulae 

M q q = (H® + + # (2) > (46) 

H<® is the sum of the binding energy £ a and the current masses m s = 150MeV, m c = 
1350MeV, mj = 4700MeV, (we take all light flavors to be massless). H^ 1 ' stands for the 
induced instanton mass for the light AM q ~ 420MeV and heavy AMq ~ 86MeV quarks. 
f/"( 2 ) provides for the extra Coulomb binding energy discussed above the 1/rriQ hyperfme 
splitting within the multiplets. Our results are summarized in Table 1. Overall, the results 



are in reasonable agreement with experiment and the constituent quark model [15]. 



B. Baryons 

Heavy baryons may be analyzed in similar fashion using the induced interactions dis- 
cussed at the end of section 3. First, we note that the analog of the one-gluon exchange in 
our case, is the induced two-body interaction ([32|), that is 

and similarly for baryons with two heavy quarks. This interaction is expected to be overall 
attractive, thus binding. We note that it scales like N c in baryonic configurations. Instantons 
in heavy baryons, induce also a three-body interaction (|38|), that is 
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x {^.[i + Iataj-^.^a?-^] 

+ \*q- [A? + ^ ■ i^ a h - ■ °i *5 ■ (A 6 A B )i + (« <- j) ] } (48) 

This interaction scales as N®. Although sub leading in our previous book-keeping arguments, 
we will keep it in our N c = 3 arguments. 

For baryons, the trial wavefunctions, will be chosen in the form 

Hx,v) = (^-f 2 e~ ai ^ 2) (49) 

where x — 73(^1 — ^2) and fj = + r 2 — 2r*g) are standard Jacobi coordinates. Here 

we choose r qq q = 1 fm for the size of the heavy-light baryons. The size of the heavier 
configurations will be set to r q QQ ~ 0.86 fm for Qqq and tqqq ~ 0.7 fm for QQQ. The 
confining energy is about S a ~ 500 MeV for all of them. 

The addition of one more quark in the baryonic configurations, brings about in the three- 
body contribution to the energy an additional overall factor of lZ q for a light quark, and IZq 
for a heavy quark, in comparison to (f47|). Specifically, 

The three-body contribution is repulsive, whereas the two-body contribution is attractive. 
Thus, there is a subtle interplay between two- and three-body interactions in the determi- 
nation of the overall energy of the heavy-light baryonic systems. For a baryon size r = 1 fm, 
TZ q ~ —0.75 and IZq ~ —0.06. We note that the results are very sensitive to the size of the 
hadron. Indeed, for r = 0.9 fm, lZ q ~ —1, and for r = 0.4 fm, TZq ~ — 1. In other words, 
two- and three-body contributions, become comparable in strength. The size r is fixed by 
the choice of the potential ( ^7J ) and is independent from the character of the induced in- 
teraction in our discussion. Here, we chose to work with r q Q = 0.6 fm for the heavy-light 
mesons, and r qq Q ~ 1 fm, for the heavy baryons. A smaller size, say r qq Q ~ 0.5 fm would 
not fit the spectrum. Of course, other choices may also be possible, with other choices of 
the potential in ((40D . 
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We present the spectra for heavy baryons in Table 2. Our results are in a reasonable 
agreement with experiment and other models. Major uncertainty comes from the 1/N C 
corrections (non-planar graphs) and the approximation for the ground-state wavef unctions. 
As we have seen, the first order corrections were up to 25% of the leading order. Since we are 
making two expansions, AQ CD /mg and l/N c , one would also expect substantial contributions 
coming from the first order corrections in 1/N C for physical N c = 3. 

C. Exotics 

The rationale of constructing two- and three-body interactions using instanton induced 
effects, can be extended to multi-quark configurations. For example, for QQqq configu- 
rations, the integration over group (cf. ( |A-9| ) ) leads to the four-body interaction of the 
form 

/AM ? \ 2 /AM Q \ 2 

(— 1 + 7° —1+7° _ _ 

x (Q — 2 — Q Q — 2 — ^ (^M^ + ^QrVl ) 

+ ^ Q^^A a Q Q^p^Q (detq L A a q R + detq R X a q L ) 

+ \ Q^^A a Q Qi±^A a Q (detq^q* + detq fi q L ) ) (50) 

The overall sign is consistent with the naive expectation, that the ra-body interaction follows 
from the (n + l)-body interaction by contracting a light quark line, resulting in an overall 
minus sign (quark condensate). 

We recall that for each extra light quark, the penalty factor in the energy is lZ q . Starting 
with r = 1 fm for three quark states (whether heavy or light), we find that the radius 
r shrinks to 0.93 fm for one additional light quark (four-quark state), to 0.88 fm for an 
extra one (five-quark state), and to 0.84 fm for still another one (six-quark state). For 
r = 1 fm, lZ q ~ —0.72, while for r = 0.9 fm, 1Z q — — 1. It follows that the three-body 
interaction (repulsive) will tend to overcome the binding energy provided by the two-body 
interaction (assumed attractive) in the multiquark configurations of the type (qqqq), {Qqqq), 
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(qq qqq), (Qq qqq) and (qqq qqq). In this respect, we agree with the conclusions of |L6| that 



the H-dibaryon viewed as a six-light-quark state (qqq qqq), will be unbound by the three 
body-forces induced by instantons 

Adding a heavy quark brings about a penalty factor TZq in the energy. This factor is 
0.06 for r = 1 fm, and 1 for r = 0.4 fm. Using the harmonic potential with two light quarks 
and four heavy quarks, yield r = 0.8 fm and TZq = 0.1. Thus, for heavy six-quark states 
the three-body interaction is 10% of the two-body interaction, hence small. In this respect, 
if we were to think about the H-dibaryon cLS cL SIX - heavy-light-quark state (Qqq Qqq) will 
not be unbound by the three-body interaction. Similarly, the four-body-interaction will 
be expected to be about 1% of the two-body, about the same order of magnitude as the 
hyperfine splitting discussed above. Therefore we conclude, that the multi-body effects are 
only important for multi-quark states near threshold. 

We have run specific calculations for exotics containing two and three heavy quarks 
of the type (QQ qq) and (QQ Qqq) . With our choice of parameters, we have found that 
these configurations were stable against strong decays through QQ qq — > Qq + Qq and 
QQ QQQ ~^ QQ + QQQ or QQ QQQ ~> QQ + QQq- In both cases the binding energy was found 
to be of the order of 10 MeV, in agreement with other models WJ\ and 



V. CONCLUSIONS 

We have presented a general framework for discussing the effects of a dilute and random 
instanton gas, on the multiquark configurations involving heavy and light flavors. Our 
instanton-induced interactions obey chiral and heavy quark symmetry to leading order in 
the bare heavy quark mass. Recoil and spin effects were explicitly worked out and found to 
be small on single quarks. 

We have used the instanton-induced interactions to analyze the spectra of heavy-light 



In Ref. [16] the radius r = 0.5 fm was used in comparison to r = 0.84 fm in our case. 
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mesons and heavy-light baryons. The results are in overall agreement with the constituent 
quark model results as well as soliton calculations. The role of the three-body force in 
multiquark states was also discussed. If the strangeness is viewed as a light degree of 
freedom, then the H-dibaryon may be unbound by three-body effects. If on the other hand, 
strangeness is viewed as a heavy degree of freedom, then the faith of the H-dibaryon is 
controlled by the strength of the two-body forces. 

We have analyzed the role of multi-body induced instanton interactions on multi-quark 
states and found them to be very sensitive to the size of the states in light systems. The 
size of the system is fixed by long range confining forces, and thus outside the scope of the 
instanton-based models. In heavy systems, the two-body interaction is dominant whatever 
the size of the systems considered. 

Clearly the present analysis could be extended in several directions. First, the derivation 
was based on the planar approximation, and that could be lifted as subdominant effects 
may be considered. Second, the spectrum calculations may be refined, by considering more 
realistic potentials and trial wavef unctions. A more thorough analysis of the correc- 
tions could be carried along the lines we have discussed using Bethe-Salpeter construction. 
In this respect, it would be interesting to re-investigate directly the present effects on the 
various correlation functions. Finally, one could use it to calculate magnetic moments and 
other static characteristics of the hadrons. 
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APPENDIX A: U - INTEGRATION 



Here we present a graphical shortcut to the derivation of the higher order interactions. 
Our task is to average over the string of color matrices 

fduf[ ut:u^ = [dU[ Uftf{ T (A-l) 
J k=l J 

with invariant measure / dU = 1. Averaging over color is equivalent to finding all projections 
onto the singlets of the group, i.e., P(3 ® 3) n — > 1. For small n — 1, 2, 3 the answer can be 
obtained by a direct calculation using the method developed by Creutz [|T^] for averaging 
over links in lattice gauge theory. For higher n the color integration using this approach 
becomes very involved and leads to cumbersome expressions involving sums of n! 2 strings 
of 2n Kronecker delta's. On the other side, since 3 <S> 3 = 1 © 8, the problem reduces in 
practice to finding all projections of the product of n octets (adjoint representations) onto 
the singlet, i.e., P(®8) n — > 1 for SU(N) with N = 3. The number of distinct projections 
A n grows with n like 

n -i 

A n = n\Y J {-l) k V] (A-2) 

In order to avoid explicit presentation of the indices, we use the diagrammatic technique, 
originally proposed by Cvitanovic |EU . 



Fundamental graph consists of links and vertices based on the following identification 
(see Fig. 1): 

V?U*i - + [At [At' (A-3) 

where A a (a = 1, . . . , iV 2 — 1) are the color Gell-Mann matrices, normalized as TrA J A J = 2 <5 lJ . 
To each projector and symbol we assign a weight as follows 

Pi = —, P 8 = r (A-4) 

1 iV 8 4(iV 2 -l) V ; 

For n = 1 one has 
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J dUUfrfi = P(l © 8) = P x = ±% 5\. (A-5) 

since averaging over one octet leads to zero. The indices on the r.h.s. follow from the 
graphical representation, n = 2 makes use of Pg, i.e., projector of two octets onto singlet 
(Fig. 2) 



J dU[U?U ]3 h } 2 = P(l <g> 1 + 8 ® 8) = Pi ■ Pi + P 8 (8 ® 8 

$ 2 + 4(jV 2 1 - i) M£ MS (A-6) 



Contributions to the spin (ij) and color (ab) parts are totally identical, and to simplify the 
notation, we retain only the color matrices in our formulae and graphs. 

For n = 3 the only new ingredient is the projection of the product of three octets onto 
singlets. Fig. 3 shows how to project n — 3. Explicitly fijM, 



P(®8) 3 = P[8 <g> (1 + 8 a + 8 S + 10 + 10 + 27)] = P 8 [8 ® 8 a + 8 <g> 8 S 
= P%Pd d(color • spin) + PgPf /(color • spin) 

' A ^ k (\\\i\ k 3 r b • d(spin) 



4(iV 2 - 1) 2(iV 2 - 4) 
+ ^ P k m^)l ■ /(spin) (A-7) 

where <i(spin) = d abc {X1\ b l X'^)l and we have used values of the symmetric and antisymmetric 
octet projectors 

1 N 
P f = w> Pd= 2(N* - 4) (A " 8) 

Generalization of the above procedure for arbitrary n is possible due to the classification 
of all invariant tensors for 577(3) by Dittner |2"I|| . All of them could be constructed from 
the combinations of f abd and d abc , standard structure constants defined by commutator 
[A a , A b ]_ = 2if abc \ c and anticommutator [A a , \ b ] + = 2d abc \ c + ^5 ab , respectively. Assigning 
a weight for each symbol: for f abc «-> Pf, and for d abc «-> Pj, we can derive the result for 
arbitrary n from the graphical representation. 
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For n = 4 one can either have two P$ projections or start with P§ and then contract 
using four combinations dd, fd, df, and // 



J dU[U?UVf = ^[6 4 } a b M J dU[U?UVf + sym. perm's 



(spin)] • (AiA J 2 A*Ai)£ { a(at : ^. a . at : 5^5 kl + 5 lk 8> 1 + S U S^ 



1 1 



4(iV 2 - l)4(iV 2 - 1) 



M i xikx.il i xilxik 



+ 4(iV 2 -l) PdPd 



d ijm d klm + tfkm d jlm + Jim ^fcmj + f. (A-9) 

Our last example is n = 5, Fig. 5. Symbolically, it reads 



P(®8) E 



P 8 • P 8 P d + P 8 ■ P d P d P d +{d^f) 



x (color • spin) (A-10) 



Any other n can be considered in a similar fashion, and since the explicit expression becomes 
lengthy, the use of the diagrams helps greatly in analyzing the properties of the derived 
formula. 

One should keep in mind that not all of possible (/, d) combinations are linearly inde- 
pendent and using appropriate relations between them ||21|| , one can reduce their number 
to the number ( |A-2| ). Our expression have an advantage of being explicitly symmetric. We 



have tested the procedure in various ways. One way is to take all possible contractions. For 
example, as a consistency check, one may contract a pair of indices, e.g., {b\a 2 ) (iij 2 ) and 
see if it is reduced to the n — 1 case 

6%6% J dU[U?uV} n = N J dU[U?uV b r l (A-ll) 

The same procedure can be used to check formulae for any n. Another useful check is 
the analysis of the powers of N. The leading behaviour must scale like 1/N n for any n 
as iV — > oo. The leading contribution in this limit has a form of the determinant build of 
Kronecker deltas, reproducing, in the instanton model case, 't Hooft's result for an arbitrary 
number of light flavors. 
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Figure Captions 

Figure. 1. Basic relation. 
Figure. 2. n = 2 contraction. 
Figure. 3. n = 3 contraction. 
Figure. 4. n = 4 contraction. 
Figure. 5. n = 5 contraction. 
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Table Captions 

Table I. Mesonic spectrum. 
Table II. Baryonic spectrum. 
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FIGURES 
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FIG. 1. Basic relation 
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FIG. 2. n = 2 contraction 
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FIG. 3. n = 3 contraction 
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FIG. 4. n = 4 contraction 




FIG. 5. n = 5 contraction 
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TABLES 



TABLE I. Mesonic spectrum 



I{J) F Exp. value Prediction QM |15| 

D° i(0~) 1869 1881 1800 -r- 1860 

D* i(l-) 2010 1965 1930 -=- 1990 

Df 0(0- ) 1969 2031 1975 

Df 0(1") 2110 2115 2061 
£° i(O-) 5278 5231 



TABLE II. Baryonic spectrum 





I 


J 


Exp. value 


Prediction 


QM H 


Ref. 


17| 







1 

2 


(2285) 


2376 


2200 


2170 




1 


1 
2 


(2453) 


2502 


2360 


2421 


— C 


i 

2 


1 
2 


2468 


2652 


2420 


2421 







1 
2 


(2704) 


2802 


2680 


2645 



? 3558 3550 3510 

? 3708 3730 3698 

? 4808 4810 4784 
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